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Glutaredoxins (Grxs) are small oxidoreductases particularly specialized in the reduction 
of protein-glutathione adducts. Compared to other eukaryotic organisms, higher plants 
present an increased diversity of Grxs which are organized into four classes. This work 
presents a thorough comparative analysis of the biochemical and catalytic properties 
of dithiol class I Grxs from poplar, namely GrxCI, GrxC2, GrxC3, and GrxC4. By 
evaluating the in vitro oxidoreductase activity of wild type and cysteine mutated variants 
and by determining their dithiol-disulfide redox potentials, pKg values of the catalytic 
cysteine, redox state changes in response to oxidative treatments, two subgroups can be 
distinguished. In accordance with their probable quite recent duplication, GrxCI and GrxC2 
are less efficient catalysts for the reduction of dehydroascorbate and hydroxyethyldisulfide 
compared to GrxC3 and GrxC4, and they can form covalent dimers owing to the presence 
of an additional C-terminal cysteine (Cysc). Interestingly, the second active site cysteine 
(CysB) influences the reactivity of the catalytic cysteine (CysA) in GrxCI and GrxC2 
as already observed with GrxC5 (restricted to A. tlialiana), but not in GrxC3 and C4. 
However, all proteins can form an intramolecular disulfide between the two active site 
cysteines (CysA-Cyse) which could represent either a protective mechanism considering 
that this second cysteine is dispensable for deglutathionylation reaction or a true catalytic 
intermediate occurring during the reduction of particular disulfide substrates or in specific 
conditions or compartments where glutathione levels are insufficient to support Grx 
regeneration. Overall, in addition to their different sub-cellular localization and expression 
pattern, the duplication and maintenance along evolution of several class I Grxs in 
higher plants can be explained by the existence of differential biochemical and catalytic 
properties. 
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INTRODUCTION 

Glutaredoxins (Grxs) are oxidoreductases sharing a conserved 
3D structure with members of the thioredoxin (Trx) superfam- 
ily. They are present in most living organisms except in some 
bacterial and archaeal phyla (Rouhier et al, 2008; Alves et al., 
2009). The primary function of Grxs was long thought to be the 
reduction of disulfide bonds and more particularly those formed 
between reduced glutathione (GSH) and a protein cysteinyl 
residue, a process known as glutathionylation, but as explained 
below, specific Grx members could rather serve as iron-sulfur (Fe- 
S) cluster transfer proteins (Rouhier, 2010). Glutathionylation 
is a post-translational modification which potentially fulfills 
several functions. It can constitute an intermediate of some 
reaction mechanisms, but it can also modulate protein func- 
tion and serve as a signaling mechanism or protect cysteine 
residues from irreversible oxidation in plants (Zaffagnini et al., 
2012b) and in animals (Dalle-Donne et al, 2009). Although a 
few additional proteins, as some specific thioredoxins, can cat- 
alyze deglutathionylation reactions, Grxs are likely the major 



deglutathionylation system, at least in plants (Bedhomme et al., 
2012; Chibani et al., 2012). They are usually regenerated via a 
NADPH/glutathione reductase (GR)/GSH system, but a few Grxs 
can be reduced by ferredoxin- or NADPH-dependent thiore- 
doxin reductases (Johansson et al., 2004; Zaffagnini et al., 2008). 
Depending on the Grx isoforms and target proteins or substrates 
considered, Grxs can use different catalytic mechanisms, called 
monothiol and dithiol mechanisms (Rouhier et al., 2008). The 
monothiol mechanism is used for the reduction of glutathiony- 
lated proteins and requires a priori only the catalytic cysteine (the 
first or more N-terminal of the two active site cysteines, referred 
to as CySA). It performs a nucleophilic attack on the protein- 
glutathione adduct, the Grx becoming glutathionylated. This 
oxidized Grx is regenerated by reduction via a GSH molecule. The 
dithiol mechanism requires the catalytic cysteine but also a recy- 
cling cysteine which could be either the second active site cysteine 
(Cyss) or an additional extra active site cysteine (Cysc)- If the 
target protein is glutathionylated, the first step is similar to the 
monothiol mechanism, but the glutathione-adduct formed on 
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Grx is solved by one of these recycling cysteines instead of GSH. 
The resulting disulfide is then either reduced by two molecules 
of GSH or by a thioredoxin reductase. If the target protein has 
an intra or inter-molecular disulfide, the mechanism is compa- 
rable to the one used by Trxs. The catalytic cysteine forms a 
transient mixed disulfide with one of the cysteines of the tar- 
get protein, which is then resolved by the recycling cysteine. As 
above, final regeneration occurs via two molecules of GSH or via 
a thioredoxin reductase. 

Based on the active site sequence, Grxs were initially cat- 
egorized into two classes, a dithiol (CP[Y/F]C motif) and a 
monothiol (CGFS motif) class (Rodriguez-Manzaneque et al., 
1999) that were subsequently renamed classes I and II, respec- 
tively (Couturier et al., 2009a). Analyses of Grx distribution and 
evolution in archaea, bacteria, and eukaryotes suggested that Grx 
domains of classes I and II have evolved through duplication 
and divergence from one initial gene present in the last com- 
mon ancestor of all organisms (Alves et al, 2009). Two additional 
classes are found in land plants (Couturier et al., 2009a). The Grx 
isoforms belonging to class III are specific to land plants and those 
from class IV are also present in few algae and animals. There 
is no biochemical or functional information on the latter Grxs 
yet. From genetic analyses, it appears that class III Grx isoforms 
also named CC-type Grxs or ROXY play a role in plant develop- 
ment and in pathogen defense mechanisms, reviewed in (Meyer 
et al., 2012). These functions are related to their ability to inter- 
act with TGA transcription factors, likely regulating their redox 
states. However, the biochemical properties of these Grxs have 
not been studied in detail yet, because of the difficulty to pro- 
duce soluble recombinant proteins (Couturier et al., 2010). Only 
one study demonstrated that Arabidopsis thaliana GrxSlO is able 
to regenerate in vitro the mitochondrial type IIF peroxiredoxin 
(Finkemeier et al., 2005). The Grxs from class II are particular 
considering their involvement in the maturation of Fe-S clus- 
ters and in the regulation of iron homeostasis (Rouhier et al, 
2010), which is likely related to their capacity to bind labile Fe- 
S clusters and to transfer them to target proteins (Bandyopadhyay 
et al., 2008). On the other hand, though many studies performed 
with class II Grxs from various sources indicated that they have 
no or very poor efficiency for the reduction of the traditionally 
tested glutathionylated substrates using a GR/GSH regeneration 
system, at least one study indicated that they have a reduc- 
tase activity. Indeed, Chlamydomonas reinhardtii Grx3 is able to 
reduce in vitro a glutathionylated A4-glyceraldehyde-3-phosphate 
dehydrogenase (Zaffagnini et al, 2008). Interestingly, it can be 
regenerated by the ferredoxin-thioredoxin reductase but not by 
the GSH/GR system. The existence of alternative regeneration 
systems may explain why the characterization of such activities 
for class II Grxs has been retarded. Overall, the demonstrated 
capacity of several class I Grxs to provide electrons to a battery 
of enzymes, as peroxiredoxins and methionine sulfoxide reduc- 
tases, enzymes that go through a catalytic cycle which involves 
glutathionylation of their catalytic cysteine (Rouhier et al., 2001; 
Vieira Dos Santos et al., 2007; Tarrago et al, 2009), suggests that 
they are the central and favorite deglutathionylating agents in 
cells. 

In plants, class I regroups six different glutaredoxin members 
called GrxCl, C2, C3, C4, SI 2, and C5, the latter being restricted 



to Brassicaceae (Couturier et al, 2011). The cytoplasmic GrxCl 
and the plastidial GrxC5 are able to bind an [Fe2S2] cluster, at 
least in vitro into recombinant proteins (Rouhier et al., 2007; 
Couturier et al., 2011). The biochemical and structural charac- 
terization of the two plastidial isoforms, GrxS12 and GrxC5, con- 
firmed that the apoforms of these Grxs reduce glutathionylated 
proteins with catalytic properties quite similar to other plant and 
non-plant class I Grxs (Couturier et al., 2009b, 2011). Moreover, 
poplar GrxS12 was shown to be sensitive and temporarily inacti- 
vated after treatments with oxidizing compounds such as hydro- 
gen peroxide (H2O2), nitrosoglutathione (GSNO), and oxidized 
glutathione (GSSG) (Zaffagnini et al, 2012a). Considering the 
very electronegative redox potential of the glutathione adduct 
formed on the catalytic CySA (— 315mV at pH 7.0), the reduc- 
tion of GrxS12 necessitates a quite low redox potential of the 
GSH/GSSG couple in the chloroplast. Thus, GrxS12 may act as 
a stress-related redox regulator allowing glutathione to play a 
signaling role in some oxidizing conditions by maintaining the 
glutathionylation of its target proteins over the duration of the 
stress period. Several class I Grxs, C. reinhardtii Grxl and Grx2 
(Gao et al, 2010), A. thaliana GrxCl, C2 (Riondet et al, 2012) 
and GrxC5 (Couturier et al., 2011) and poplar GrxCl, GrxC4, 
and GrxS12 (Rouhier et al., 2002, 2007; Couturier et al, 2009b; 
Zaffagnini et al., 2012a) have been partially or extensively char- 
acterized at the biochemical level, but some members have not 
yet been studied and there is no thorough and comparative study 
allowing to understand if the existence of several members is only 
related to their sub-cellular localization or to their expression pat- 
tern or whether some biochemical properties can partially explain 
their diversity. 

From the biochemical analyses performed here with the four 
dithiol class I Grxs existing in poplar, GrxCl to GrxC4 (the 
fifth member, GrxS12, has a monothiol active site and was char- 
acterized previously), two sub-groups with distinct properties 
emerged. In addition, this work provides new insights on the 
existence of different oxidation forms involving the active site 
cysteines (CySA and Cyss) or an additional C-terminal cysteine 
(Cysc) and on their potential role for the reaction mechanism. 

MATERIALS AND METHODS 
SITE-DIRECTED MUTAGENESIS 

The construction of pET3d expression plasmids containing 
poplar GrxCl, C2, C3, and C4 was described previously (Rouhier 
et al, 2002, 2007). All cysteine residues of GrxC2, C3 were indi- 
vidually substituted into serines using two complementary muta- 
genic primers (listed in Table SI). The corresponding mutated 
proteins are called GrxC2 C23S, C26S, or C80S, GrxC3 C37S or 
C40S. The recombinant proteins GrxCl C31S, C34S, C88S, and 
GrxC4 C27S and C30S were prepared previously (Rouhier et al., 
2002, 2007). Two additional mutants, in which only the catalytic 
cysteine is remaining, GrxCl C34/88S and GrxC2 C26/80S, have 
also been cloned. 

HETEROLOGOUS EXPRESSION AND PURIFICATION OF RECOMBINANT 
PROTEINS IN £ coli 

For protein expression, cultures of the E. coli BL21(DE3) strain, 
containing the pSBET plasmid and co-transformed with the 
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different recombinant pET3d plasmids, were successively ampli- 
fied up to 2.4 L in LB medium at 37°C supplemented with 
50 (JLg/mL of ampicUlin and kanamycin. Induction of pro- 
tein expression was performed at exponential phase by adding 
100 (jlM isopropyl P-D-thiogalactopyranoside for 4h at 37°C. 
After centrifugation (20 min at 4400 g), the cell pellets were resus- 
pended in about 20 mL of TE NaCl buffer (30 mM Tris-HCl 
pH 8.0, 1 mM EDTA, 200 mM NaCl) and eventually conserved 
at — 20°C. Cell lysis was performed by sonication and the solu- 
ble and insoluble fractions were separated by centrifugation for 
45 min at 27,000 g. 

Purification of proteins was carried out in three steps. The sol- 
uble fraction was first precipitated by ammonium sulfate from 0 
to 40% and then to 80% of the saturation. The 40-80% ammo- 
nium sulfate-precipitated fraction was subjected to gel filtration 
chromatography (ACA44 gel) equilibrated with TE NaCl buffer. 
Afl:er dialysis against TE (30 mM Tris-HCl pH 8.0, 1 mM EDTA) 
buffer and concentration, the interesting fractions were loaded 
to a DEAE (diethylaminoethyl) sepharose column equilibrated in 
TE buffer. GrxCl WT, GrxC2 WT and their corresponding cys- 
teinic mutants passed through the column, whereas GrxC3 WT, 
GrxC4 WT and their respective cysteinic mutants were retained 
and eluted using a linear 0-0.4 M NaCl gradient. The purest frac- 
tions, as judged by SDS-PAGE analysis, were pooled and dialyzed 
against TE buffer by ultrafiltration in Amicon cells equipped with 
a YMIO membrane. Finally, the fractions were concentrated and 
stored at — 20°C until further use. Protein purity was checked by 
SDS-PAGE, and protein concentrations were determined spec- 
trophotometrically using the corresponding molar extinction 
coefficients at 280 nm of 10,095 M^^ cm^^ for GrxCl WT and 
its monocysteinic mutants, 9970 M^^ cm^' for GrxCl C34/88S, 
8605 M^^ cm^' for GrxC2 WT and its monocysteinic mutants, 
8480 M^i cm^i for GrxC2 C26/80S, 7575 M^^ cm^^ for GrxC3 
WT, 7450 M^i cm"' for GrxC3 C37S and C40S, 4595 M^^ cm"' 
for GrxC4 WT and 4470 M^^ cm^^ for GrxC4 C27S and C30S. 

ENZYMATIC ACTIVITIES 

The thioltransferase activities [2-hydroxyethyl-disulfide (HED) 
and dehydroascorbate (DHA) assays] were measured as described 
previously (Couturier et al., 201 1) using as purified non-reduced 
proteins. Briefly, measurements were performed at 25°C in 
steady-state conditions by following NADPH oxidation at 340 nm 
in the presence of a Grx reducing system which is composed 
of NADPH, GR, and GSH. The reaction was started by adding 
Grx after a 3-min pre-incubation time and Grx activity was cor- 
rected by subtracting the spontaneous reduction rate observed 
in the absence of Grx. Because of the spontaneous reaction of 
GSH with the substrates, we can generally not work at saturating 
concentrations of all substrates using these assays. Enzyme and 
substrate concentrations used are indicated in the legend of fig- 
ures and tables. The activity was expressed as nmol of NADPH 
oxidized/nmol of Grx/s using a molar extinction coefficient of 
6220 M"^ cm"^ at 340 nm for NADPH. Three independent 
experiments were performed at each substrate concentration, and 
kcat and apparent K,„ values were calculated by non-linear regres- 
sion using the Michaelis-Menten equation using the program 
GraphPad Prism 4. 



PREPARATION AND ANALYSIS OF OXIDIZED Grxs 

Around lOmg of protein was reduced using 40 mM DTT in 
500 |xl of 30 mM Tris-HCl pH 8.0 buffer for 2 h at 25°C followed 
by desalting on G25 column pre-equilibrated with 30 mM Tris- 
HCl pH 8.0 buffer. Oxidized proteins were prepared by incubating 
500 [xM reduced poplar Grxs in ca 200 (jlL with a 10 fold excess 
GSSG or H2O2 or a 5 fold excess GSNO for 2h at 25° C before 
desalting on G25 columns. The protein oxidation state was then 
analyzed by electrospray ionization mass spectrometry analysis as 
described previously (Couturier et al., 201 1). 

For alkylation assays, reduced Grxs treated or not with oxi- 
dants were diluted to 10|xM into 50|xl of 30 mM Tris-HCl pH 
8.0 buffer and precipitated on ice for 30 min with one volume 
of 20% trichloroacetic acid (TCA). After centrifugation (10 min 
at 13,000^) and washing with 100% acetone, the pellet was 
resuspended into 10|xl of 100 mM Tris-HCl pH 8.0, 1% SDS 
containing 2 mM of methoxyl-PEG maleimide of 2 kDa (mPEG 
maleimide) which alkylates free thiol groups. The protein mixture 
was then separated on non-reducing 15% SDS-PAGE. 

For midpoint redox potential titrations, CySA glutathionylated 
forms were obtained by the above described GSNO treatment of 
the double cysteinic mutants (GrxCl C34/88S, GrxC2 C26/80S) 
and of the monocysteinic mutants (GrxC3 C40S and GrxC4 
C30S). The intramolecular disulfide bond (CysA-Cyss) contain- 
ing proteins were obtained by the above described H2O2 treat- 
ment but using GrxCl C88S, GrxC2 C80S and GrxC3 and GrxC4 
WT isoforms. 

pKa DETERMINATION OF THE CATALYTIC CYSTEINE 

The pKa measurements of N-terminal active site cysteines (CySA) 
have been performed with WT Grx isoforms following a proce- 
dure described in (Gallogly et al, 2008). Briefly, 3 [iM reduced 
Grx was incubated with or without 300 |i,M iodoacetamide (lAM) 
in 100 mM sodium citrate or phosphate buffers ranging from pH 
2.0 to 7.0. Following this pre-incubation step, Grx activity was 
determined by adding an aliquot of the pre-incubation mixture 
to the HED assay described above. The percentages of remaining 
activity at each pH were determined by comparing the activity of 
the enzyme incubated with and without JAM and an adaptation 
of the Henderson-Hasselbach equation (Gallogly et al., 2008) was 
used for pKa value determination using the program GraphPad 
Prism 4. 

MIDPOINT REDOX POTENTIAL (f m) DETERMINATION 

Oxidation-reduction titrations using the fluorescence of the 
adduct formed between protein free thiols and monobromobi- 
mane (mBBr) were carried out at ambient temperature in 200 |i,l 
of 100 mM HEPES pH 7.0 buffer containing 10|xM of oxi- 
dized proteins, either glutathionylated or with an intramolecular 
disulfide-bond, and defined mixtures of oxidized and reduced 
DTT or glutathione for more positive values to set the ambient 
potential (£;,). Total concentration of DTT or glutathione was 
2 mM. After 2 h incubation, mBBr was added at a final concen- 
tration of 2 mM and the reaction was carried out in the dark 
for 20 min. The reaction mixture was then precipitated on ice 
for 30 min with one volume of 20% TCA. After centrifugation 
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(lOmin at 13,000 g) and washing with 100% acetone, the pel- 
let was resuspended into 400 (il of 100 mM Tris-HCl pH 8.0, 
1% SDS. Fluorescence emission of the resulting solution was 
measured at 472 nm after excitation at 380 nm using a variant 
Gary Eclipse (Agilent). Values were transformed into percent- 
ages of reduced protein and fitted to the Nernst equation using 
non-linear regression for Em value determination. 

RESULTS 

KINETIC COMPARISON OF POPLAR CLASS I Grxs 

An amino acid sequence alignment of the predicted mature forms 
of all poplar class I Grxs reveals several conserved residues, par- 
ticularly those around the active sites and those involved in glu- 
tathione binding (Figure 1). However, GrxGl and GrxG2 present 
GGYG and GPFG active site motifs, respectively, different from 
the GPYG active site found in GrxG3 and GrxG4. Interestingly, 
contrary to GrxG3 and GrxG4, GrxGl and GrxG2 isoforms pos- 
sess an additional cysteine (Gysc) located in their G-terminal part 
in a IGGGD motif, as in human Grxl, E. coli Grx3 and the two 
other class I Grxs found in plants, GrxC5 and GrxS12 (Figure 1). 



With the aim of comparing the redox and kinetic properties of all 
dithiol Grxs (GrxC5 is not found in poplar) and understanding 
why these duplicated genes have been retained during evolu- 
tion, we have performed a thorough biochemical analysis of the 
corresponding recombinant proteins. 

The catalytic properties of each isoform were first determined 
by measuring their reductase activity using the classical HED and 
DHA assays, which measure the reduction of glutathionylated- 
P-mercaptoethanol or of dehydroascorbate, respectively. All pro- 
teins display a better catalytic efficiency (kcatlKm) in the HED 
test with a difference of two orders of magnitude (Table 1). In 
both assays, the kinetic parameters determined for each protein 
follow the same tendency. GrxG3 and G4 have quite compa- 
rable turnover numbers, with k^at values about 10 fold higher 
than those of GrxGl and G2. In terms of catalytic efficiencies, 
the dichotomy between GrxGl/G2 and GrxC3/G4 is clear but 
it is often below this factor of 10 because GrxG3 and G4 have 
higher apparent _K",„ values for both substrates. Overall, GrxGl 
always exhibited the lowest catalytic efficiency and GrxG4 the 
highest. Goncerning the apparent affinity for GSH, although there 
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FIGURE 1 I Amino acid alignment of class I glutaredoxins. Tlie 
alignment was done using ClustalW from sequences of the 
recombinant proteins that are deprived of putative targeting 
sequences. The strictly conserved amino acids are depicted in white 



on black. Other conservative amino acid changes are indicated in black 
on gray. Residues involved in GSH binding based on the 3D 
structures of poplar GrxC4 and GrxS12 and Arabidospis GrxCS are 
marked by an asterisk. 
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Table 1 | Kinetic parameters of poplar dithiol class I glutaredoxins. 







PtGrxCI 


PtGrxC2 


PtGrxC3 


PtGrxC4 


PtGrxS12' 


AtGrxC5= 


p-Mt-oo 


Km (mM) 


U. lU ± U.U 1 


U.Ub ± U.U 1 


U.oU ± U.Ub 


U.oz ± U.U4 


U.o 1 ± U.U4 


U.zU ± U.UZ 




1^ /r-— 1 \ 


7.9 ± 0.2 


11 .7 ± 0.1 


77.0 ± 4.4 


on o 1 /I n 

yu.o ± 4.y 


TO ^ inn 

zo. 1 ± u.y 


1 o 1 1 n no 
1 .Z 1 ± U.Uo 




WCmlxloBM-i.s-i) 


0.79 


1.95 


1.54 


2.82 


0.745 


0.06 


DHA 


Km (mM) 


0.55 ±0.06 


0.37 ±0.03 


1.14±0.13 


0.60 ±0.09 


0.38±0.19 


0.21 ±0.03 






0.6 ±0.1 


1.3 ±0.1 


6.5 ±0.4 


8.2 ±0.6 


1.7 ±0.2 


0.23 ±0.01 






1.09 


3.51 


5.70 


13.67 


4.6 


1.1 


GSH 


Km (mM) 


2.76 ±0.26 


3.87 ±0.53 


3.60 ±0.42 


1.96 ±0.28 


4.0 ±0.5 


nd 




kcatis-'') 


13.8 ±0.6 


30.1 ±2.3 


131 .3 ±8.3 


102.9 ±6.5 


92.8 ±5.3 


nd 




kcat/Kmix W^M-''.s-'') 


5.00 


7.78 


36.47 


52.50 


23.20 


nd 



The apparent Km values for fi-ME-SG (glutathionylated (i-mercaptoethanol) were determined using a IHED concentration range of 0. 1-1 mM in presence of 2mM 
GSH and 100 nM enzyme for GrxCl, 50nM for GrxC2, WnM for GrxC3 and GrxC4. The apparent Km vaiues for DHA were determined using a DHA concentration 
range of 0.1-1 mM in presence of2mM GSH and 2^1. M enzyme for GrxC 1, Iv-Mfor GrxC2, 200 nM for GrxC3, and 150 nM for GrxC4. The apparent Km values for 
GSH were determined using a GSH concentration range of 0.25-5mM in the presence of 0.7mM HED and lOOnM enzyme for GrxCI, SOnM for GrxC2, lOnM 
for GrxC3 and GrxC4. The feat and Km values were calculated by non-linear regression using the Michaelis-Menten equation. Values are the means ± SD of three 
replicates. "Data extracted from the following references Couturier et al. (2009b, 2011), Zaffagnini et al. (2012a). nd, not determined. 



are differences in the K„, values for GSH between the four iso- 
forms, ranging from ~2.0 mM for GrxC4 to ~3.9 mM for GrxC2, 
the GrxCl/C2 and GrxC3/C4 dichotomy is not visible using this 
parameter. 

Then, the involvement for the catalytic mechanism of the dif- 
ferent cysteine residues present in each isoform was investigated 
by comparing the activity of WT proteins to their corresponding 
cysteinic mutants using similar assays at fixed substrate concen- 
trations (Figure 2). As expected, all monocysteinic mutants for 
the catalytic cysteine CySA (GrxCl C31S, GrxC2 C23S, GrxC3 
C37S, and GrxC4 C27S) were inactive, confirming that it is indis- 
pensable for the reductase activity. On the contrary, mutating 
Cysc, in positions 88 and 80 in GrxCl and GrxC2, respectively, 
did not affect catalytic activity. Finally, the mutation of CySB has a 
different impact depending on the isoforms. The protein variants 
for this cysteine in GrxCl and C2 (GrxCl C34S and GrxC2 C26S) 
were 5 fold more efficient and 2-3 fold more efficient than the 
corresponding WT proteins in the HED and DHA assays, respec- 
tively. On the other hand, the mutation of this cysteine in GrxC3 
and GrxC4 has almost no impact, since, only a slight decrease was 
observed in the HED assay for the GrxC3 C40S variant. 

OXIDATIVE MODIFICATIONS OF CYSTEINE RESIDUES 

Several observations prompted us to analyze the sensitivity of 
all dithiol class I Grxs to oxidizing treatments and the resulting 
changes in their redox state. First, the mutation of CysB increased 
the catalytic efficiency of GrxCl and C2 but not of GrxC3 and 
C4 suggesting possible difference in the reaction mechanisms. 
Second, although Cysc does not influence their catalytic activ- 
ity, the migration of the purified proteins in non-reducing SDS 
PAGE gels indicated the formation of covalent linkages for GrxCl 
and C2 that are not visible for GrxC3 and C4 (data not shown). 

For this purpose, reduced proteins have been treated for 2 h 
with a 10 fold excess of GSSG or H2O2, or a 5 fold excess of 
GSNO. The resulting redox state was analyzed by mass spectrom- 
etry and by protein migration on non-reducing SDS PAGE with 
or without thiol alkylation by 2kDa mPEG maleimide, which 



allows distinguishing reduced and oxidized forms of Grxs. Since 
reduced mPEG-alkylated proteins and covalent dimers of GrxCl 
and GrxC2 migrate roughly at the same position, controls with- 
out alkylation have been included (Figure 3). With non-alkylated 
proteins, we confirmed that an H2O2 treatment, but not the 
GSSG and GSNO treatments, promoted the formation of cova- 
lent dimers of GrxCl and GrxC2 (Figures 3A,B) but not of 
GrxC3 and GrxC4 (Figures 3C,D). Only a small part of GrxC2 
dimer is visible after these 2 h GSSG and GSNO treatments. From 
the procedure that included an alkylation step, we can further 
conclude that all treatments led essentially to a complete oxi- 
dation of all thiol groups present on all Grxs since almost no 
shift of migration was observed. A minor form with one free 
thiol group is visible for GrxC2 after a GSSG treatment. It is 
even more discrete upon GSNO treatment. For GrxC3 and C4, 
these oxidized monomeric forms are interpreted as containing 
intramolecular disulfides between active site cysteines (CySA- 
Cysfi) (Figures 3C,D). For GrxCl and C2, these two cysteines 
should also form this intramolecular disulfide, and Cysc is either 
modified by a glutathione or a nitroso (in the case of GSNO) 
adduct as confirmed by mass spectrometry analyses (see below) 
or it is involved in an intermolecular disulfide upon H2O2 treat- 
ment. From similar experiments performed with the GrxCl C88S 
and GrxC2 C80S mutants (Figures 3E,F), the absence of covalent 
dimers indicated that Cysc is indeed involved in dimerization. 
Moreover, in the context of these variants, after GSSG and GSNO 
treatments, the proteins mainly contain an active site disulfide 
bond whereas a small part of the proteins has one remaining free 
thiol, likely corresponding to glutathionylated intermediates. 

In order to identify the type of protein modification formed 
by these oxidative treatments, all proteins have been analyzed by 
mass spectrometry (Table 2 and Figures S1-S6). First, all reduced 
proteins gave a single species with a mass decrease of ca 131 Da 
compared to the calculated theoretical molecular masses. It does 
undoubtedly correspond to the cleavage of the first methionine 
as expected from the presence of an alanine in second position. 
The results obtained with GrxC3 and GrxC4 are more easily 
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FIGURE 2 I Reductase activity of GrxC1, C2, C3, and C4 and of 
cysteine-to-serine variants. (A) HED reduction was measured in the 
presence of 2 mM GSH and 0.7 mM HED using either 100 nM GrxCI , 
20 nM for GrxCI C34S, 50 nM GrxC2, 10 nM forGrxC2 C26S, 10 nM GrxC3, 
or 10 nM GrxC4. (B) DHA reductase activity was measured in the presence 
of 2 mM GSH and 1 mM DHA using either 2 ixM GrxCI , 1 |iM for GrxCI 
C34S, 1 |xM GrxC2, 250 nM GrxC2 C26S, 200 nM GrxCS, or 150 nM GrxC4. 
The data are expressed as nmol NADPH oxidized per nmol enzymes per s 
and represent means ± SD of three separate experiments. 



interpreted since these proteins contain only two cysteines. In 
the presence of H2O2, and in accordance with the absence of 
detected free thiols in the previous experiments (Figures 3C,D), 
GrxC3 and GrxC4 are exclusively detected as a monomeric oxi- 
dized form with a decrease of ~2 Da confirming the formation 
of a CysA-CysB intramolecular disulfide bond. The same forms 
were primarily obtained upon GSSG and GSNO treatments but a 
minor monomeric form with a mass increment of 305 Da corre- 
sponding to one glutathione adduct was also detected (Table 2). 
It was not detected on gels presumably because of the presence of 
minute amounts. Nevertheless, this suggests that the glutathiony- 
lated form represents an intermediate for the formation of the 
intramolecular disulfide. 



The results obtained for GrxCI and GrxC2 are slightly more 

complex. In response to an H2O2 treatment, two different species 
are formed for both proteins. In accordance with the results from 
alkylation experiments, a major form corresponds to a disulfide- 
bridged dimer having also an intramolecular CysA-CysB disulfide 
bond. The detected monomeric forms, which have no free thi- 
ols or a very small fraction as observed for GrxC2, exhibited a 
mass increment of 29.5 Da for GrxCI and 32.2 Da for GrxC2, 
suggesting the presence of two oxygen atoms. Considering that 
this overoxidation is not observed for GrxCI C88S and GrxC2 
C80S variants which only form monomeric proteins with an 
intramolecular disulfide (as indicated by the mass decrease of 
2 Da), the results obtained for the WT proteins are best inter- 
preted as the formation of an intramolecular disulfide between 
the active site cysteines (CySA-CysB) and of a sulfinic acid on 
Cysc- In the presence of GSSG or GSNO, GrxCI and GrxC2 exist 
only under monomeric forms but with variable oxidation states. 
Except for the GSSG-treated GrxCI, there is a minor form with 
a mass decrease of about 2 Da consistent with the formation of 
an intramolecular disulfide. Another minor form, found in all 
samples except the GSNO-treated GrxC2, exhibited a mass incre- 
ment of 610.5-611 Da, which corresponded to the presence of 
two glutathione adducts. However, the major species formed by 
these two treatments is a monomer with a mass increment of 
about 303.5 Da which is consistent with the presence of one glu- 
tathione adduct on Cysc in addition to the CysA-CysB disulfide 
bond. Accordingly, the major species detected in protein variants 
mutated for Cysc (Grx CI C88S and GrxC2 C80S) is the one 
with an intramolecular disulfide (decrease of 2 Da upon GSSG or 
GSNO treatments). The glutathionylated forms detected by mass 
spectrometry were visible on gels (see Figures 3E,F). 

The in vitro observation of intermolecular disulfide bonds 
in GrxCI and C2, involving Cysc, raised the question of their 
possible reduction. Thus, we tested the ability of DTT but also 
of a physiological GSH reducing system composed of NADPH, 
GR, and GSH to reduce this intermolecular disulfide, the forma- 
tion of which was initiated by a 2h H202treatment. As shown 
in Figure 4, both DTT and GSH are able to reduce GrxCI and 
GrxC2 dimers to a similar extent. Whether minute amounts of 
reduced monomeric Grx initially present or formed can catalyt- 
icaUy increase the reduction of dimeric Grxs in the presence of 
GSH cannot be excluded. Overall, it indicates that the formation 
of these disulfide-bridged forms of GrxCI and GrxC2 isoforms 
is reversible, notably by GSH, a relevant cellular physiological 
reductant. 

REDOX PROPERTIES OF CLASS I GLUTAREDOXINS 

As in other thiol-dependent oxidoreductases, the catalytic effi- 
ciency of class I Grxs should be at least partially governed by 
thermodynamic parameters as the redox midpoint potentials of 
dithiol-disulfide couples and by the reactivity of the catalytic cys- 
teine which is dependent on its p_fCa. Hence, we determined the 
CysA pJCfl of the four Grx isoforms using a method relying on 
iodoacetamide which is an alkylating reagent reacting with thio- 
lates but not thiols. The pH-dependent inactivation of WT Grx 
by iodoacetamide was followed by comparing the activity of WT 
proteins in the HED assay, after pre-incubation of the protein in 
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FIGURE 3 I Sensitivity of poplar Grxs to oxidative treatments. 

Pre-reduced Grxs were subjected to 2 h oxidative treatments with H2O2, 
GSSG, and GSNO as described in the "materials and methods" section. Prior 
to separation on SDS-PAGE gels in non-reducing conditions, GrxCl WT (A), 
GrxC2 WT (B), GrxC3 WT (C), GrxC4 WT (D), GrxCl C88S (E), and GrxC2 
C80S (F) were alkylated or not with 2 kDa mPEG maleimide. The shift 



observed following the alkylation of the free thiol groups is larger than 
expected (ca. 4 kDa instead of 2 kDa), but this has been already observed 
previously (Chibani et al.. 2011 ). The numbers on the right correspond to the 
number of alkylated thiols that remained free after reducing or oxidizing 
treatment of Grxs. The stars in panels (A) and (B) indicate the position of 
covalent dimer. 



different buffers ranging from pH 2.0 to 7.0 in the presence or not 
of iodoacetamide. From these titration curves, we obtained pKa 
values of 5.3 ± 0.1 for GrxCl, 5.0 ± 0.1 for GrxC2, 4.6 ± 0.1 for 
GrxC3, and 4.6 ±0.1 for GrxC4 (Figure 5). These results indicate 
first that at physiological pH, the thiolate ion will be the dominant 
species for all Grxs. Although the differences are not huge, it also 
confirmed the existence of two subgroups, GrxCl and C2 having 
higher p^C^compared to GrxC3 and C4. 

Two disulfide bond forms may be catalytically relevant for all 
these Grxs, the one with a glutathionylated CySA and the one 
with a CySA-CysB intramolecular disulfide bond. For measur- 
ing the midpoint redox potentials (£„) of each of these disul- 
fide/dithiol couples, we used different protein variants and two 
distinct oxidative treatments in order to discriminate between 
these two oxidized forms. To obtain proteins containing an 
intramolecular disulfide bond, reduced Grxs possessing only the 
two active site cysteines (GrxCl C88S, GrxC2 C80S, GrxC3 WT, 
and GrxC4 WT) have been treated with H2O2 before desalting 
on G25 columns. On the other hand, glutathionylated proteins 
have been obtained by treating reduced proteins in which only 
CysA was remaining (GrxCl C34/88S, GrxC2 C26/80S, GrxC3 
C40S, and GrxC4 C30S) with GSNO before desalting. The val- 
ues determined at pH 7.0 for the intramolecular disulfide range 



from —233 mV for GrxC3 to —263 mV to GrxCl (Figure 6). On 
the other hand, only little variations were visible when measur- 
ing E,„ values for the glutathione adduct formed on the catalytic 
cysteine. They varied from — 242mV for GrxC3 to — 254mV 
for GrxCl. Overall, the intramolecular active site disulfide bond 
and glutathionylated forms of GrxC2 and GrxC4 present simi- 
lar Em values. A 10 mV difference was measured between both 
forms for GrxCl and C3, the more electronegative redox poten- 
tial being the active site disulfide in GrxCl and the glutathione 
adduct in GrxC3. However, the differences are small and should 
not explain the observed differences in catalytic efficiencies. In 
reducing cellular environments, all these oxidation forms should 
be compatible with glutathione reduction. Considering their pos- 
sible targeting to secretory pathways, this may not be the case for 
Grx3 and C4 if they are in the endoplasmic reticulum or in the 
apoplasm. 

DISCUSSION 

DUPLICATION AND NEO-FUNCTIONALIZATION INTO CLASS I 
GLUTAREDOXINS 

Among all living organisms, higher plants have the highest num- 
ber of genes coding for Grxs. In addition to specific isoforms 
such as those belonging to the class III, land plants also possess 
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FIGURE 4 1 GSH-dependent reduction of the intermolecular disulfide in 
poplar GrxC1 and GrxC2. The ability of GSH to reduce poplar 
disulfide-bridged dimer GrxCl and disulfide-bridged dimer GrxC2 generated 
by a 2 h H2O2 treatment was evaluated in the presence of 1 mM DTT and 
1 mM GSH coupled to NADPH and GR. The lanes Ox represent the initial 
oxidized proteins. Proteins were separated on SDS-PAGE gels under 
non-reducing conditions. 
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FIGURE 5 I p^a determination of the catalytic cysteine of poplar 
GrxCI, C2, C3, and C4. Reduced GrxCl (A), GrxC2 (B), GrxC3 (C), and 
GrxC4 (D) WT proteins were incubated in different buffers ranging from pH 
2.0 to 7.0 in presence or not of iodoacetamide (lAM) prior to measurement 
of their activity by using the HED assay. The percentages of remaining 
activity at each pH were determined by comparing the activity of the 
enzyme incubated with and without lAM and an adaptation of the 
Henderson-Hasselbach equation was used for pKg value determination. 
The data are represented as mean ± SD of three separate experiments. 



an expanded class I, usually with 5-6 isoforms (Couturier et al., 
2009a). From previous phylogenetic analyses, it was quite clear 
that, there are three Grx couples, GrxCl/C2, GrxC3/C4, and 
GrxC5/S12 that probably evolved by more or less recent dupli- 
cation (Couturier et al., 2009a). Interestingly, GrxCl and C2 
proteins are cytoplasmic proteins (Riondet et al., 2012), GrxC3 



E„s.s=-250±1mV 
= -247 ± 1 mV 




E„s-s=-2'»6±1 mV 
EmssG=-250±1mV 



I » s-s o sss I 

FIGURE 6 I Determination of redox potential of oxidized poplar GrxCl, 

C2, C3, and C4. Representative curve of oxidation-reduction titration of 
CysA-CysB intramolecular disulfide forms (black circles) and glutathionylated 
CySft (white circles) forms of poplar GrxCl (A), GrxC2 (B), GrxC3 (C), and 
GrxC4 (D) at pH 7.0. The titration was carried out using a total DTT or GSH 
concentration of 2 mM for 2 h and resulting free thiol groups were labeled 
by mBBr. Values in the text are the means ± SD of three replicates. 



and C4 have N-terminal extensions predicted to target the pro- 
teins to secretory pathways but the final destination is not known, 
and GrxC5 and S12 are chloroplastic (Couturier et al, 2009b, 
2011). The case of GrxCl and C2 is interesting. Whereas GrxC2 
prototypes are found in all land plants (mosses, lycophytes, 
monocots, dicots), GrxCl is specifically found in dicots. A redun- 
dancy between both isoforms is suggested in A. thaliana by the 
observation that single knock-out mutants in grxcl or grxc2 are 
aphenotypic, whereas a double mutant is lethal at an early stage 
after poUinization (Riondet et al., 2012). Considering that GrxCl, 
but not GrxC2, is able to bind an Fe-S cluster (at least in vitro), 
this would rather suggest that the essential nature of these two 
paralogs is linked to their reductase activity. These results point 
to the need of investigating at the genetic level a possible redun- 
dancy among the GrxC3/C4 and GrxC5/12 couples. This study 
primarily focused on the biochemical properties of dithiol class 
I Grxs since a detailed structure-function analysis was previously 
achieved on the two other prototypes of class I Grxs, the poplar 
monothiol GrxS12 and the Arabidopsis dithiol GrxC5 (Couturier 
et al., 2009b, 2011). GrxC5 isoforms are only found in species 
of the Brassicaceae family. We have demonstrated in particular 
that, although being very similar in sequence due to their recent 
duplication, a single substitution in the active site sequence (from 
WCSYS in GrxS12 to WCSYC in GrxC5) explains the capacity of 
GrxC5, but not GrxS12, to bind an [Fe2-S2] cluster into a dimer. 
Among other class I Grxs, only GrxCl can also exist under a holo- 
form binding an Fe-S cluster (Rouhier et al., 2007). In this case, 
the mutational analysis showed that the presence of the glycine in 
the CGYC active site instead of a proline in other Grxs (CPFC in 
GrxC2 or CPYC in GrxC3 and C4) was sufficient to explain this 
property. Overall, this indicates that subtle sequence differences 
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affecting key residues may provide different properties to other- 
wise quite close proteins. This prompted us to further analyze the 
various glutaredoxins of this class as GrxC2 and GrxC3 have been 
only poorly studied so far. 

TWO SUBGROUPS IN DITHIOL Grxs EXHIBIT DISTINCT KINETIC 
PROPERTIES 

Considering the divergence observed in the amino acid sequences, 
the first aim of this v/oik was to perform a thorough compar- 
ative analysis of the biochemical properties of GrxCl, GrxC2, 
GrxC3, and GrxC4 from poplar. All these proteins exhibit deg- 
lutathionylation or reductase activities using classical HED and 
DHA assays with catalytic efficiencies (10^ M^' s^^ for HED 
and 10^-10^ M^' s^^ for DHA) in the same range as those pre- 
viously reported for plant and non-plant class I Grxs (Table 3). 
Two subgroups can however be distinguished, GrxCl and GrxC2 
appeared to be less efficient catalysts than GrxC3 and GrxC4. 
The difference is particularly clear when considering the turnover 
numbers (kcat) and the kcat/^^DHA and kcat/^^GSH- This is not vis- 
ible with kcatlK^-ME-SG because the differences in the turnover 
number are compensated in this case by differences in apparent 
affinity in favor of GrxCl and C2. An observation that may be 
interesting to understand these differences is that GrxCl/C2 and 
GrxC3/C4 should have quite different charge distribution at the 
protein surface that could account for substrate recognition and 
thus for the observed differences in catalytic efficiency. Indeed, 
their behavior on ion exchange chromatography indicates that 
GrxCl and GrxC2 should be more basic having a pi close or 
higher than 8, whereas GrxC3 and GrxC4 should be more acidic 
having a pi value below 8 (theoretical pi estimations are indeed 
comprised between 5 and 6). 

In order to further explain the differences between poplar class 
I Grxs, we determined the pKa value of their catalytic cysteine, 
as it was proposed that reactivity of human and Chlamydomonas 
class I Grxs is partially related to the pKa of N-terminal active 
site cysteine (Gallogly et al., 2008; Gao et al, 2010). The values 
obtained for poplar Grxs are globally close to those reported for 
other plant and non-plant class I Grxs (Table 3), including the 
value of 3.9 determined for poplar GrxS12 using the same pro- 
cedure (Zaffagnini et al, 2012a; Roos et al., 2013). Interestingly, a 
difference of ca 0.5 unit is visible between GrxCl/C2 (pKg around 
5) and GrxC3/C4 (pKa around 4.6). Such differences between 
paralogs have been described already since C. reinhardtii Grxl, 
S. cerevisiae Grx2, and human Grxl exhibit higher turnover num- 
bers often translated into higher catalytic efficiencies compared to 
C. reinhardtii Grx2, S. cerevisiae Grxl, and human Grx2, respec- 
tively (Tables). For the human and C. reinhardtii enzymes, it 
is also associated to a difference in pKa values of ca 1 unit, the 
enzymes with the lowest pKa being the most efficient. According 
to the Bronsted theory, in a thiol-disulfide exchange reaction, each 
unit decrease in the pKa value of the leaving group is responsible 
for a 4 fold increase in the second-order rate constant (Szajewski 
and Whitesides, 1980). This could fit with the conclusion that dif- 
ferences in pKa values between class I poplar Grx isoforms are 
at least partially responsible for the differences in catalytic effi- 
ciencies observed. On the other hand, only small differences were 
observed for the oxidation-reduction midpoint potentials of both 



the glutathione adduct or the intramolecular disulfide and they 
are not discriminating GrxCl/C2 vs. GrxC3/C4. The measured 
values (from —230 to —260 mV) are in the same range than those 
reported for E. coli Grxl (—233 mV) and human Grxl and 2 
(-232 and -221 mV) (Table 3) (Aslund et al., 1997; Sagemark 
et al, 2007). 

THE ROLE OF THE RECYCLING Cyss: MONOTHIOL OR DITHIOL 
REACTION MECHANISM? 

Interestingly, the differences in protein activities between GrxCl 
or C2 and GrxC3 or C4 are attenuated when CysB of GrxCl or 
C2 is removed (see Figure 2). This has been reported also for 
Arabidopsis GrxC5, yeast Grxl, human Grxl and Grx2, and pig 
Grx (Yang and Wells, 1991; Yang et al, 1998; Johansson et al, 
2004; Gallogly et al., 2008; Discola et al, 2009; Couturier et al, 
2011). In clear contrast, mutating this cysteine in yeast Grx2, E. 
coli Grxl and Grx3 led to a decrease of Grx activity (Bushweller 
et al, 1992; Nordstrand et al, 1999; Peltoniemi et al., 2006; 
Discola et al., 2009; Saaranen et al., 2009). For E. coli Grxl and by 
extension for Grxs for which mutation of this cysteine decreased 
their activity, this decrease was attributed to a change in glu- 
tathione recognition, binding, or affinity, because CysB would 
determine the glutathione specificity of the glutathionylated Grx 
reduction step (Saaranen et al, 2009). For GrxCl and GrxC2 and 
other eukaryote Grxs mentioned above in which the mutation of 
this cysteine increased their activity, the data indicates that the 
presence of CysB slows down the reaction for some reason. It may 
eventually modif)' the pKa of CySA. This has been experimentally 
confirmed for a pig Grx, for which the pKa of the catalytic Cys22 
is 3.8 in a WT protein but 4.9 and 5.9 in C25S and C25A variants 
(Yang and Wells, 1991). Another possibility is that a CySA-CysB 
intramolecular disulfide is formed during activity measurements 
either because it constitutes a regular step of the reaction cycle 
in some situations, especially for dithiol substrates or because it 
is a side reaction reflecting the competition between CySB and 
GSH for the reduction of the glutathionylated CysA (Gallogly 
et al., 2009). Then, the reduction of this CysA-CysB intramolec- 
ular disulfide intermediate would require two molecules of GSH 
vs. one molecule of GSH for the glutathionylated form and add 
two steps in the reaction mechanism. Thus, it seems unfavorable 
in terms of efficiency and energetic cost. 

We have previously analyzed possible variations in the redox 
state of the dithiol AtGrxC5 (WCSYC active site) in response to 
15 min treatments with GSSG or oxidized DTT by measuring the 
tryptophan intrinsic fluorescence of the protein, since the oxi- 
dation of the catalytic cysteine quenched the fluorescence of the 
adjacent tryptophan (Couturier et al, 2011). It was concluded 
that the dithiol AtGrxC5 cannot form a CySA-Cysn intramolec- 
ular disulfide first because only the GSSG treatment led to 
the formation of a glutathionylated protein and to fluorescence 
quenching and second because the WT protein was crystallized 
with a glutathionylated adduct. However, at that time, we did use 
neither H2O2 and GSNO treatments, nor prolonged incubation, 
nor alkylation assays. In fact, reconsidering the results of mass 
spectrometry showing, in addition to the peaks with one or two 
glutathione adducts, a peak with a decrease of 2 Da between the 
reduced protein and the GSSG-treated protein, it is very likely 
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that a small fraction of the protein contained an intramolecular 
disulfide formed between active site cysteines (Couturier et al, 
2011). Here, we have clearly observed that a prolonged oxidative 
treatment of reduced class I Grxs with H2O2, GSSG, or GSNO 
led preferentially to the formation of an intramolecular disul- 
fide bond between CysA and CySB- For the H2O2 treatment, it 
likely goes first through the formation of a sulfenic acid on CySA 
before CysB performs a nucleophilic attack on it. For the GSSG 
or GSNO treatments, it is very likely that a nucleophilic attack 
of CysA onto these oxidized glutathione species leads to the for- 
mation of a glutathione adduct on Cysyi which is then slowly 
reduced by a Cys^ nucleophilic attack leading to the formation 
of the intramolecular disulfide and to the concomitant release of 
a GSH molecule. 

It is worth mentioning here that these in vitro experiments 
have been done in particular with the aim of exploring all pos- 
sible post-translational modifications but the conditions are not 
exactly physiological first because the concentrations of the oxi- 
dants may not be realistic and more importantly because this 
is done in the absence of GSH, which is normally found at 
high amounts and is highly reduced in most cell compartments. 
The presence of GSH would certainly prevent to a large extent 
the formation of these intramolecular disulfides because it may 
preferentially reduce either the sulfenic acids or the glutathione 
adducts over the Cyss recycling cysteine. This is certainly the case 
for GrxC3 and GrxC4 since mutating CysB has no effect on the 
activity. The formation of an intramolecular disulfide in these 
proteins could thus constitute uniquely a side reaction, protect- 
ing the cysteine or inactivating the proteins only under specific 
conditions. On the contrary, if we consider that the increase of 
protein activity observed for CysB mutated GrxCl and GrxC2 
in HED and DHA assays in the presence of a fully reduced 
GSH pool reflects the formation of this intramolecular disul- 
fide, it may be physiologically relevant for GrxCl and GrxC2. 
This is interesting in several respects. First, while the current 
dogma is that the preferential substrates of Grxs are glutathiony- 
lated proteins, recent studies performed with mammalian Grx2 
indicated that they can reduce dithiol substrates, e.g., proteins 
with intra- or intermolecular disulfide (Hanschmann et al, 2010; 
Schutte et al., 2013). In this case, the Grxs should employ a reac- 
tion mechanism similar to the one used by Trxs involving both 
CySA and CySB. Second, based on the observation that some 
Grxs can be reduced by thioredoxin reductases (Johansson et al., 
2004; Zaffagnini et al, 2008; Couturier et al, 2013), another 
possibility could be that in specific physiological situations or 
sub-cellular compartments where glutathione is depleted, absent 
or oxidized, the formation of an intramolecular disulfide may 
be determinant either to favor the use of an alternative reducing 
system as the one constituted by NADPH and thioredoxin reduc- 
tase or to simply constitute a protective mechanism, preventing 
the irreversible overoxidation of this residue and thus protein 
inactivation. 

In conclusion, the role of CySB in dithiol Grxs is yet uncertain 
and may depend on the Grx isoform as this residue influences dif- 
ferentially the protein activity. The fact that the formation of an 
intramolecular disulfide bond between CySA and CysB is possible 
in all cases lets open the possibility to use both monothiol and 



dithiol catalytic mechanisms, contrary to Grxs possessing only 
one cysteine residue. 

THE ROLE OF THE SEMI-CONSERVED Cysc IN GrxCl AND GrxC2 
ISOFORMS 

The Cysc cysteine residue, present in the IGGCD motif (at posi- 
tion 88 and 80 in GrxCl and C2, respectively), is found in several 
other class I Grxs, including plant GrxCS and S12 as well as in 
human Grxl and E. coli Grx3. Moreover, it is also present in 
many class II Grxs. It was shown for instance that a CySA-Cysc 
intramolecular disulfide bond can be formed in S. cerevisiae Grx5 
and C. reinhardtii Grx3 (Tamarit et al, 2003; Zaffagnini et al, 
2008), suggesting that Cysc can serve as a recycling cysteine in 
the absence of CysB- Reduced glutathione does not seem to be 
an efficient reductant for this intramolecular disulfide bond and 
in the case of CrGrx3, it is reduced by a ferredoxin-thioredoxin 
reductase system (Zaffagnini et al., 2008). Although GrxS12 has a 
monothiol WCSYS active site sequence, a CysA-Cysc intramolec- 
ular disulfide has not been observed whatever the oxidizing 
conditions tested. Moreover, in the crystal structure of a glu- 
tathionylated GrxS12, the sulfur atom of Cys87 is 9.6A away from 
the Cys27-SG adduct (Couturier et al., 2009b). Hence, this cys- 
teine may not serve as a recycling cysteine in class I Grxs and does 
not a priori participate to the reaction mechanism. Consistently, 
its mutation in GrxCl and C2 did not influence protein activity, as 
previously observed for poplar GrxS12, Arabidopsis GrxCS, or E. 
coli Grx3 (Nordstrand et al, 1999; Couturier et al, 2009b, 2011). 

However, considering its position very close to the active 
site residues and to residues involved in GSH binding an effect 
on activity could have been observed. Indeed, it is positioned 
between the so-called GG kink which is important in determin- 
ing the backbone geometry of the following amino acids, and 
a Thr or an Asp which forms hydrogen bonds with GSH (see 
Figure 1) (Couturier et al, 2009b; Li et al, 2010). In the GrxS12 
structure, the backbone amino group of Cysc also directly forms 
a hydrogen bond with the glutamyl group of GSH (Couturier 
et al., 2009b), but the substitution of the cysteine into ser- 
ine would not disrupt this interaction, preventing to definitely 
rule out a GSH stabilizing effect of this cysteine. On the other 
hand, it is probable that modifying the redox state of this cys- 
teine might have more pronounced effects influencing for exam- 
ple GSH binding. It has been shown for example that, under 
non-reducing conditions, E. coli Grx3 exists under an expected 
monomeric form but also forms disulfide-bridged dimer or mul- 
timers via Cysc (Aslund et al., 1994, 1997). In addition, previously 
reported oxidative treatments of human Grxl, which contains 
three extra active site cysteines including Cysc, led to the iden- 
tification of several possible post-translational modifications, i.e., 
intramolecular disulfide, disulfide-bonded dimers and oligomers, 
glutathione adducts, nitrosylation, some of them inhibiting its 
activity (Hashemy et al., 2007). Depending on the conditions, 
Cysc was found to be either nitrosylated or glutathionylated 
or involved into a disulfide. Similarly, we have observed redox 
changes of GrxCl and C2 upon treatment with oxidants that 
could further provide clues about the possible function of this cys- 
teine, if any. In the presence of H2O2, this cysteine was found to be 
either overoxidized or involved in disulfide-bridged homodimers. 



Frontiers in Plant Science | Plant Physiology 



December 2013 | Volume 4 | Article 518 | 12 



Couturier et al. 



Biochemical characterization of poplar dithiol glutaredoxins 



Since this dimerization only occurs in presence of H2O2, the for- 
mation of an intermediate sulfenic acid is likely required before 
intermolecular disulfide bond formation which is consistent with 
the presence of a small part of overoxidized proteins. Interestingly, 
GSH can efficiently reduce these disulfide-bridged dimers, mak- 
ing this reaction physiologically reversible. In the presence of 
GSSG or GSNO, Cysc is also prone to oxidative modification, but 
in the form of a glutathione adduct. Although we cannot test the 
impact of these modifications on Grx activity (we cannot selec- 
tively modify Cysc without also oxidizing CySAand CySfi), we can 
speculate that the different reversible oxidation forms of Cysc 
observed here might represent regulatory or signaling intermedi- 
ates transiently modifying the functions of these Grxs by affecting 
for example the nature or the time-course reduction of the target 
proteins. 
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Figure S1 | Electrospray ionization mass spectrometry analysis of PtGrxCI 

WT. Spectra of whole protein were determined for reduced protein before 
(A) and after treatment with GSSG (B), GSNO (C), or H2O2 (D) as 
described in the Methods section. 

Figure S2 | Electrospray ionization mass spectrometry analysis of PtGrxC2 

WT. Spectra of whole protein were determined for reduced protein before 
(A) and after treatment with GSSG (B), GSNO (C), or H2O2 (D) as 
described in the Methods section. 

Figure S3 | Electrospray ionization mass spectrometry analysis of PtGrxCS 

WT. Spectra of whole protein were determined for reduced protein before 
(A) and after treatment with GSSG (B), GSNO (C), or H2O2 (D) as 
described in the Methods section. 

Figure S4 | Electrospray ionization mass spectrometry analysis of PtGrxC4 

WT. Spectra of whole protein were determined for reduced protein before 
(A) and after treatment with GSSG (B), GSNO (C), or H2O2 (D) as 
described in the Methods section. 

Figure S5 | Electrospray ionization mass spectrometry analysis of PtGrxCI 
C88S. Spectra of whole protein were determined for reduced protein 
before (A) and after treatment with GSSG (B), GSNO (C), or H2O2 (D) as 
described in the Methods section. 

Figure S6 | Electrospray ionization mass spectrometry analysis of PtGrxC2 
C80S. Spectra of whole protein were determined for reduced protein 
before (A) and after treatment with GSSG (B), GSNO (C), or H2O2 (D) as 
described in the Methods section. 

Table S1 | Primers used for site-directed mutagenesis. 
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